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Martensitic transformation plays a pivotal role in the microstructural evolution and plasticity of many 
engineering materials. However, so far the underlying atomic processes that accomplish the displacive 
transformation have been obscured by the difficulty in directly observing key microstructural signatures on 
atomic scale. To resolve this long-standing problem, here we examine an AISI 304 austenitic stainless steel 
that has a strain/microstructure-gradient induced by surface mechanical attrition, which allowed us to 
capture in one sample all the key interphase regions generated during the y(fcc) £(hcp) a'(bcc) 
transition, a prototypical case of deformation induced martensitic transformation (DIMT). High-resolution 
transmission electron microscopy (HRTEM) observations confirm the crucial role of partial dislocations, 
and reveal tell-tale features including the lattice rotation of the a' martensite inclusion, the transition lattices 
at the £/a' interfaces that cater the shears, and the excess reverse shear-shuffling induced y necks in the £ 
martensite plates. These direct observations verify for the first time the 50-year-old 
Bogers-Burgers-Olson-Cohen (BBOC) model, and enrich our understanding of DIMT mechanisms. Our 
findings have implications for improved microstructural control in metals and alloys. 



The structural transformation involved in fast cooling of iron and steel was first studied by Adolf Martens in 
the late 19*^ century. In honor of his contributions, a category of displacive crystalline structural phase 
transitions has since been termed "martensitic transformation" \ Martensitic transformations are diffusion- 
less and characterized by a collective movement of atoms across distances that are typically smaller than one 
nearest-neighbor spacing. The resulting structural transition is often accompanied by pronounced changes in 
mechanical (e.g. strength, toughness, and shape memory^"^), electromagnetic (e.g. magneto- and mechano- and 
electro -caloric^), and transport (such as conductivity^) properties of the materials. In biological systems, some 
virus species use pressure-induced martensitic transformation to infect bacteria cells^'^. Martensitic transforma- 
tions of silicates in the Earth's mantle are thought to play a vital role in plate tectonics and earthquakes 

The importance of martensitic transformations in nature and engineering has led to tremendous efforts to 
explore their microscopic mechanisms. Almost enjoying the equal popularity as the classic Kurdjumov & Sachs^^ 
and Nishiyama^^ shear models, in 1964, Rogers and Burgers^^ proposed a hard-sphere model for the martensitic 
transformation in steels from a face-centered cubic fcc-y austenite phase to a body-centered cubic bcc-oc' 
martensite phase. Olson and Cohen^^"^^ further developed the Rogers and Rurgers model by emphasizing the 
role of plastic deformation, proposing that the austenite y phase transits to a hexagonal close-packed hcp-s 
martensite phase first and then to the bcc-a' martensite phase. Hereafter, this model will be referred to as the 
Rogers-Rurgers-Olson-Cohen model (or the RROC model). (Supplementary Fig. 1) However, half a century later, 
the shear movements during the y ^ s ^ a' transformation hypothesized in the RROC model have yet to be 
confirmed by experiments at the atomic level. What has been known so far is that certain crystallographic 
orientation relationships (ORs) exist between the parent and martensite phases, and the resulting crystals exhibit 
fme-scale inhomogeneities involving shuffling, slip, twinning, and stacking faults^^"^^. A major obstacle in cla- 
rif)^ing experimentally the mechanism of strain-induced martensitic transformation lies in the fact that several 
processes are expected to occur in parallel in both of the parent and product phases, generating a complex pattern 
of interface dislocations and diffuse tilt grain boundaries (GRs) that are difficult to correlate and analyze. 
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Figure 1 | Orientation of bcc-a' martensites. (a), TEM images taken at the top surface layer of 304 stainless steels after SMAT. High density of bcc-a' 
formed at the intersection region of two hcp-s plates, (b-f), TEM images and the corresponding live EFT spots patterns taken along the zone axis 
[TlO]y/ /[1120]e/ /[lTl]a'. Different from the a' martensite in (b), that in (c) has grown into two hcp-s plates. The green and white dash lines in (c) 
indicate the a'main body with Pitsch orientation relationships (OR), with the parts grown into s transforming into two Kurdjumov- Sachs (K-S) ORs, 
with respect to the fcc-y. (g), Schematic summarizing the diffraction spots taken from adjoining phases in the interfacial transition regions between y, s 
and a'. The black arrows indicate the three ORs with details on parallel planes and the red arrows indicate the 5.26° misorientation angle between the 
Pitsch OR and either of the two K-S ORs. Scale bars, 200 nm (a), 5 nm (b), 10 nm (c). 



Our experiments reported below have captured the early stages of 
martensite formation. This was accomplished by inducing the y ^ s 

ol' DIMT in AISI 304 y austenitic stainless steel, via surface mech- 
anical attrition treatment (SMAT)'^ (See Methods). The SMAT 
caused a plastic strain field, which was extremely severe at the sur- 
face, with a gradient towards to the center of the sample matrix. 
Therefore, various features of the y ^ s ^ a' transformation evolu- 
tion can be captured at the different locations, along the cross sec- 
tional direction in a SMATed sample. The freshly- formed interphase 
interfaces, including their atomic arrangements with diffuse tilt GBs 
and dislocations, can then be investigated in detail under HRTEM in 
correlation with neighboring phases. The high-resolution observa- 
tions are also compared with atomistic Molecular Dynamics (MD) 
simulations (See Methods). The unprecedented level of details we 
have uncovered experimentally provides systematic evidence that 
verifies the BBOC model for the first time, and enriches our under- 
standing of the DIMT mechanism. 

Results and discussion 

Specifically, in a SMATed sample we observe high densities of dis- 
location grids in the interior of the sample matrix, the formation of 
hcp-s plates in a layer at a depth of —90 jim from the surface, and 
finally the bcc-a' martensite inclusions emerging at the intersection 
region of two s plates (Fig. la) (Supplementary Fig. 2). Once 
nucleated, the a' martensite inclusions grow towards to the 



neighboring 8 plates (Fig. Ib-c). We first analyze the ORs between 
the y, s and a' phases, which have been determined using the 
fast Fourier transformation (EFT) spots (Fig. Id-f) and atomic- 
resolution images (next section) from adjoining phases in the inter- 
facial transition regions and schematically summarized in Fig. Ig. 
The a' inclusion is found to initially exhibit the Pitsch OR with 
respect to y, when it nucleates in the intersection region between 
the two s. The Pitsch OR in our observations is different from the 
widely accepted Kurdjumov- Sachs (K-S) ORs, which require the 
closest packing plane {1 1 l}y unrotated between y, s and a' crystalline 
structures. However, when a' inclusion grows along the two 8 plates, 
the orientations are changed from the Pitsch OR to two K-S ORs in 
the grown parts. The misorientation between the Pitsch OR and 
either of the two K-S ORs is 5.26°. 

Our atomic-resolution images discern how the above is accomp- 
lished at the atomic level for the a' inclusion. Most of the HRTEM 
images in this paper are Fourier filtered to enhance the visibility 
(better contrast), at little expense of the fidelity (see typical unfiltered 
images in Supplementary Fig. 6). Fig. 2 shows the HRTEM 
images and MD simulation results on the interfaces between the 
y and a' phases and between the two ORs (all views along 
[Il0]y/ /[lIl]a'direction). We first assume that all dislocation lines 
are perpendicular to the image planes, i.e., the dislocation lines 
are parallel to the incident electron beam direction of 
[Il0]y//[1120]£//[lll]a'. The HRTEM image in Fig. 2a and the 
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Figure 2 | Interfacial structures for two orientation relationships. All projections are along[110]y//[lll]a'. The dislocations are indexed by the "±" 
symbol and its Burgers vector is depicted using Burgers circuit (dash line circuits), (a-d), HRTEM images and MD simulations of the y/a' interface with 

Pitsch OR (a,c) and K-S2 OR (b,d), respectively. (a,c), The Pitsch OR interface with (220)y//(121)a'. Two sets of dislocations b^'^ = ^ [101] on (lll)y 

Uf - - 2 

planes and b^^^ = y [101] on (1 1 l)y planes, distribute with an average distance of 8 (lOl)a' and 9 (002)y planes to release the lattice misfit strain between 

them. (b,d), The K-S2 OR interface with (lll)y//(110)a', requiring the rotated (lll)y to rotate back from the Pitsch OR. To achieve this, another 
two sets of dislocations, with Burgers vectors of -b^^^ and -b^^\ respectively, are generated to form b^^^ and -b^^^ and b^^^ and -b^^^ dislocation dipoles 
(dash-lined "_L" symbol dipoles with opposite sign). The number of -b^^^ dislocations is almost the same as the number of b^^^ dislocations, with only one 
un-coupled b^^^ dislocation, thereby making the (llT)y unrotated to evolve to the (llO)a'. The single b^^^ on every 19*^ (111)7 plane releases the lattice 
misfit strain of 5% between (lll)y and (llO)a'. (e-g), HRTEM image taken at an ol'/ol' interface, showing a 5.26° rotation inside one a' inclusion, from 
the Pitsch OR to K-S OR along the (21 1) a' interface plane, (e). Two sets of dislocations walls h^^} = y [in]on ( 110)a' planes and h^"^ = y [1 1 1] on (lOl)a' 
planes, form the diffuse small-angle tilt grain boundary. More dislocations are located in the left region of zone I, where the average dislocation spacing D 

of 12.5 A meets the requirement in the OR change, (f, g). The enlarged image in two pink-color framed rectangular zones I and II in (e), respectively. No 
dislocations are found in these regions and the lattice is continuously bending-mode deformed, leading to the diffuse interface region bent obviously and 
the (2ll)a'plane rotated 5.26°. All scale bars, 2 nm. 



MD simulation result in Fig. 2c show that the interface between 
(220)y and (121) a' with the Pitsch OR is semi- coherent and 
diffuse, in the presence of two sets of misfit dislocations deter- 
mined as b^^) = —[101] on (lll)y planes, and b^'^ = —[101] on 
(lll)y planes. These Burgers vectors have the components 

^looil=^'{?i|=^{om] = 1-8 A along the [001]y and - Z'jJio, = ^'{^iol 
= 1.276 A along the [110]y. To completely release the lattice misfit 
strain of/= 10.94% between (002)y and (IlO)a', the spacing, of 
the misfit dislocation array determined by 1^ = ^[ooi] //should be 16.5 
A, which is almost equivalent to the sum of the interplanar spacing of 



either 8 (1 01) a' or 9 (002)y planes, as shown in Fig. 2a and c. 
Furthermore, the observed dislocation number ratio for the two sets is 
closed to 1:1, which greatly reduces the long-range stress field of dis- 
locations because the b^^^ and b^^^ misfit dislocations have opposite sign 
for their Burgers vector component along the [110]y direction. 

Fig. 2b and Fig. 2d show a y/a' diffuse and semi-coherent interface 
with the K-S2 OR of (lll)y//(110)a'. When the K-S OR is evolved 
from the Pitsch OR, the rotated (lll)y must rotate back, which 
requires no component of Burgers vector along the [lll]y and 
[110] a 'directions, as shown in Fig. 2d. To achieve this, another 
two sets of dislocations are generated, which possess the Burgers 
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Figure 3 | BBOC "3T/8-T/3" shearing and evolution routes, (a), Schematic diagram illustrating the lattice evolution from the original fcc-y lattice, (i) to 
hcp-s lattices, (ii) to y' 3T/8 and y' t/3 transition lattices, and finally (iii) to bcc-a' lattice. The zone axis is[TlO]y/ /[1120]e/ /[lTl]a'. Excess T/8 and T/6 
dislocations reversely shear back to two s lattices, resulting in y' 3T/8 and y' t/3 transition lattices, respectively, satisfying the BBOC "3T/8-T/3" shearing. 
The intersecting angle of 70.52° between original (lll)y and (lll)y, is reduced to 64.08° in the y'jr/s transition lattice and 64.77° in the y'7/3 
transition lattice. The distance between neighboring atomic rows on the (lll)y and two {llljy' are 2.21 A, and the corresponding rows on {Ollja' are 
expanded to 2.34 A. (b-c), HRTEM images taken at the Si/a' (b) and 82/0^' (c) interface regions. The reduced angle of 6.44° and 5.75°, and the y' 3T/8 and 
y' T/3 transition lattices in the enlarge insets (rectangle in (b) and (c), respectively) are indicated, accomplishing the "3T/8-T/3" shearing from the two s. 
All scale bars, 1 nm. 



vectors of -b^^^ and -h^^\ respectively. The number of -b^^^ dislocations 
is almost the same as the number of b^^^ dislocations, as shown in 
Fig. 2b with only one un-coupled b^^^ dislocation, thereby recovering 
the lattice rotation caused by the b^^^ dislocation set. Fig. 2b illustrates 
that the number of -b^^^ dislocations is less than the number of b^^^ 
dislocations, leaving one b^^^ dislocation on every 19*^ (lll)y plane. 
The coupled b^^^ and -b^^^ and b^^^ and -b^^^ are dislocation dipoles, 
which do not cause any lattice rotation, thereby making the (lll)y 
unrotated to evolve to the (1 10)a' . The single b^^^ dislocation on every 
19* (lll)y plane is to release the misfit strain of 5% between the 
(lll)y and (llO)a' (See more discussions in Supplementary Fig. 3-4). 

How can the Pitsch OR and K-S OR co-survive in a a' particle 
embedded inside a single y crystal? The a' lattice must rotate itself 
inside the a' particle to fulfill the OR change. Fig. 2e shows that inside 
one a' particle, the (2ll)a' planes, or equivalently, the (Oll)a' 
planes, rotate 5.26° about [lIl]a^ leading to the a' orientation 
change from the Pitsch OR region to the K-S OR region. A small- 
angle tilt GB is formed between the Pitsch and K-S OR regions. 
Fig. 2e indicates that the tilt GB is diffuse and composed of two sets 

of dislocation walls, b^!^ = — [ill] on the (llO)a' planes and 
b^!^^ = — [111] on the (lOl)a' planes. The two sets of dislocations 



have the same edge Burgers vector component of ^[211] = 
^[211] ^ ^[2?I] " ^-^^ ^^^^g [21 1] a' direction, with which the 
dislocation spacing D = 12.8 A is anticipated from D = ^[211] Z^- 
Dislocations are not equi-distantly distributed in the two sets of 
dislocation walls. There are no dislocations at all in two pink-color 
framed rectangular zones I and II in the HRTEM image of Fig. 2e. 
More dislocations are located in the left region of zone I, where the 
average dislocation spacing D is 12.5 A, very close to the anticipated 
value of 12.8 A, meaning that the crystal lattice in the left region has 
rotated 5.26° already and meet the requirement in the OR change. 
Zooming into zones I and II we observe a severe continuous bending- 
mode deformation of the lattice, as shown in Fig. 2f and 2g, which 
makes the (2ll)a' plane rotate 5.26°. No dislocations are found in this 
region and the lattice is continuously deformed, leading to the 
obviously bent diffuse interface region. The bending-mode deforma- 
tion is elastic and interacting with the long-range dislocation stress 
field and other unknown defects under the mechanical constraints of 
the surroundings of the diffuse tilt boundary^^'^^. 

Having worked out the origin of the orientation change of the a' 
martensite and the accommodating dislocations, we now move on to 
discuss shear movements involved in the phase transition. There are 
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Figure 4 | Effect of reverse- shearing dislocations on £. (a), HRTEM images taken at the 82/0^' interface region, i.e., T/3 shearing part side. Once 
six T/6 excess shearing dislocations combine to form a reverse Shockley partial dislocation, the 82 would be transformed back to y lattice. White lines 
outline 8 atomic (1 1 l)y planes (4 atomic planes at each side), exactly one-third of the 24 (0002)8 planes, that are reverse -transformed in the neck zones at 
the y/82 coherent interface to penetrate into the 82 plate, (b). The enlarged images taken at the neck zone in (a). The stacking order of the atomic planes 
changes in the 82 plate, which is caused by the reverse Shockley partial dislocation. Scale bars, 2 nm (a), 1 nm (b). 
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corresponding to the T/3 and 3T/8 transformation dislocations, 
respectively (T donates a Shockley partial dislocation for fee twin 
shear) ^^"^^'^^ (Supplementary Fig. 1). We show next the equivalent 
of 3T/8 and T/3 shearing, in the case of the two intersecting s plates. 
The formation of a s plate can be regarded, on average, as resulting 
from T/2 transformation dislocations on every {lll}y plane. Excess 
T/6 and T/8 dislocations are therefore produced. These excess dis- 
locations can spread back to the two intersecting s plates, in a reverse 
shear-shufQing process (reverse shear of the dislocations, together 
with accommodating shuffle displacements on some planes) In 
addition to the larger volume of oc', the extremely severe strain/stress 
field caused by SMAT at the surface layer might be the driving force 
for the reverse shear- shuffling of the excess dislocations. The auste- 
nite-martensite interface might not be atomically flat, where a trans- 
ition lattice zone could exist to continuously change the lattice 
structure^^'^^'^^. Here, the reverse shear-shufQing results in two trans- 
ition lattices, which will be named y' st/s and y' t/^ transition lattices in 
the present work, with two sets of uniformly distorted {1 1 Ijy'planes^^. 
Fig. 3a schematically presents the lattice evolution from the original 
fcc-y lattice, to hcp-s lattices, to y' st/s and y' t/^ transition lattices, 
and finally to bcc-a' lattice, with the projection zone axis of 
[110]y//[1120]£//[lll]a'. The intersecting angle of 70.52° between 
the original (lll)y and (lll)y, changes to 64.08° in the y' st/s trans- 
ition lattice and 64.77° in the y' t/s transition lattice, respectively. The 
distance between neighboring atomic rows on the {lll}y, {llOlje 
and two {1 1 l}y' are all the same, 2.21 A, which must expand to 2. 34 A 
for the corresponding rows on {01 l}a'. Fig. 3b-c are, respectively, the 
HRTEM images taken at the Si/a' and 82/0^' interface regions with 
the Pitsch OR (unfiltered images in Supplementary Fig. 6). There 
exists a shuffle region between the transition lattice zone and the 
hcp-s plates. The shufQe region in the hep -82 plate is bigger than 
that in the hcp-81 plate, as the T/6 excess shearing is larger than the 
T/8 excess shearing. The angle difference of 6.44° and 5.75° (from 
70.52°) and the y' st/s and y' t/s transition lattices (see enlarged insets, 
rectangle in Fig. 3b and 3c, respectively), are thus indicative of the 
"3T/8-T/3", for shearing in the presence of the two s plates. 



We now illustrate further how the T/8 and T/6 excess dislocations 
function in the shear induced martensitic transformation. Once eight 
T/8, or six T/6, reversely spread to s and combine to form a (reverse) 
Shockley partial dislocation, y lattices or stacking faults (SFs) would be 
formed^\ Normally, SFs or extremely thin y strips are easily generated 
to accompany the hcp-81 plate, via the "T/8" action. However, the 
hcp-S2 plate is more easily transformed back to more perfect y lattices, 
because T/6 excess dislocations could combine to form more mobile 
Shockley partial dislocations. The reverse Shockley partial dislocations 
change the stacking order in the 82 plate (Supplementary Fig. 5) at the 
y/82 coherent interface, forming a y neck in the hcp-82 plate, as shown 
in Fig. 4a-b. There are 8 atomic (lll)y planes in the small y lattice 
neck zone that penetrates into the 82 plate. The number of planes is 
exactly one-third of the 24 (0002)£ planes or 24 (llO)a' planes. When 
24 (0002)£ planes convert into 24 (llO)a' planes, four Shockley partial 
dislocations must reversely slip back, which form the necks of 8 
atomic {lll}y planes. The y necks in the 82 plates are frequently 
observed at the 82/a' interface regions (Supplementary Fig. 7). 

To recapitulate, the atomic-scale observations above demonstrate 
the y/ot' Pitsch OR in the intersection region between two 8 plates, 
the y' sT/8 and y' t/s transition lattice zones at the 81/a' and 82/a' 
interfaces, and the excess dislocations involved in reverse shear. 
Taken together, the high-resolution evidence verifies the BBOC 
"3T/8-T/3" model of the strain induced y ^ 8 ^ a' transformation. 
When the yield strength of the y matrix is lower than the minimum 
stress required for martensitic transformation, plastic deformation 
will occur first in the y matrix^^ mediated by dislocations, especially 
by partial dislocations in low SF- energy metals. Depending on the 
loading rate and temperature, the partial dislocation mediated plastic 
deformation can be accomplished by the y ^ 8 phase transformation 
with T/2 dislocation slipping on every {111} plane. The subsequent 8 

a' transformation is depicted in Fig. 5 schematically. Once the 8 
plates are available in the y matrix, the formation of the a' martensite 
becomes energetically much easier. This is because the hep- 8 plates 
have already accomplished the T/2 shear displacements in the BBOC 
"3T/8-T/3" shearing model. With respect to the T/2 8 plates, the "3T/ 
8-T/3" shear displacements are reduced to the "T/8-T/6" reverse 
shearing, leading to the formation of the y' st/s and y'r/3 transition 
lattice zones. Furthermore, the excess "T/6" can cause the necks in 
the hcp-8 plates. 
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Figure 5 | Dislocation behavior in the £ ^ a' transformation. Schematic 
diagram illustrating that, with respect to the already-finished T/2 s plates, 
the BBOC "3T/8-T/3" shearing is reduced to the reverse shearing of excess 
"T/8-T/6". Reverse action of excess "T/8-T/6" transformation dislocations 
results in the SFs and y strips inside the s plates and the necks of the s plates. 

In summary, our HRTEM observations clarify the atomic arrange- 
ments and crystalline defects associated with deformation induced y 
^ s ^ a' phase transformation. The findings verify and enrich, for 
the first time, the half- century-old BBOC "3T/8-T/3" model. 
Specifically, we showed misfit dislocations at the y, s and a' interfaces 
to release misfit strains between these phases, and the diffuse tilt GBs 
that cause the change in crystal orientation of the a' inclusions. 
Furthermore, the atomic-scale images reveal the y' st/s and y'r/3 trans- 
ition lattices within the diffuse s/a' interfaces. The plastic deformation 
induced s plates reduce the "3T/8-T/3" shear to a "second step", i.e., 
the reverse action of excess "T/8-T/6" transformation dislocations, 
which result in the SFs and y strips inside the s plates and the necks 
of the s plates. The detailed understanding of the DIMT process 
elucidated here may stimulate new ideas to judiciously control the 
various ingredients/stages and the degree of transformation, and 
hence open new possibilities for the improvement of steels and alloys. 

Methods 

Sample preparation and TEM observation. The AISI 304 stainless steels used in this 
study has a composition, in wt%, of 0.04C, 0.49Si, 1.65Mn, 16.8Cr, 0.37Mo and 
balance Fe. The geometry of 304 stainless steels specimen is 70 mm X 50 mm X 
1 mm. Both surfaces of specimens were synchronously SMATed for 15 min at room 
temperature with a vibrating frequency of 20 kHz and 3 mm diameter bearing balls^*. 
The SMATed 304 stainless steels at different depth was characterized by transmission 
electron microscopy (TEM, JEOL 2010F), with operating voltage of 200 kV. 

Molecular dynamics (MD) simulations. The LAMMPS software^^ was employed in 
all MD simulations. The box size of the simulation system for the Pitsch OR was 100 
X 100 X 100 A, with bcc-a' and fcc-y phases occupying the lower, [ — 50, 0] A, and 
upper, [0, 50] A, halves, respectively, along the z direction. The lattice constants of bcc 
and fee crystals are 2.87 A and 3.602 A, respectively. When the distance between two 
atoms in fcc/bcc interface was smaller than 1.2 A, one of the atoms was removed. The 
simulation system contained 85,306 atoms. An embedded-atom method (EAM) 
potential for Fe^*^ was employed, which favors the bcc structure in the temperature 
range of 0 to 1820 K. In order to study the interfacial structure, only atoms located 
between — 3 A and 3 A of the z direction (total 4,529 atoms) were allowed to relax with 
all other atoms fixed. Periodic boundary conditions were applied along the x and y 
directions. With the NVT ensemble at temperature 1 K, energy minimization with 
conjugate gradient method was performed and followed by 1000 steps relaxation with 
the time step of 1 fs. For the K-S orientation, the initial configuration was built 
manually according the BBOC model, with 86,303 atoms in a box of 101 X 101 X 101 
A. The system was relaxed for 500,000 steps (time step of 1 fs) in NVT ensemble at 
10 K under the condition that fee atoms four layers away from the interface were 
fixed. The software of OVITO^^ and the common neighbor analysis method^" were 
employed to visualize and distinguish the atomic arrangements. 
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